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Abstract. —As more wide-range phylogenetic studies are available, the opportunity arises to 
compare DNA from these data sets to suspected introduced individuals in order to confirm species 
identification and determine their geographic origins. Two recently collected Pituophis specimens 
in Miami-Dade County, Florida, were examined using molecular analyses. Maximum likelihood 
and Bayesian inference methods place our specimens within the P. catenifer sayi I P. ruthveni 
clade. Additional morphological evidence support their identification as the Louisiana pinesnake, 
Pituophis ruthveni Stull 1929, a species indigenous to a small area in western Louisiana and eastern 
Texas and candidate for listing by the U.S. Fish and Wildlife Service. Although P. ruthveni is viewed 
as a distinct species from P. catenifer sayi based on allopatry and differences in color pattern, no 
molecular evidence was found supporting the recognition of P. ruthveni as a separate species. 
However, adding other mtDNA and nuclear DNA genes might provide needed data for distinguishing 
between these two named taxa. 
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Introduction 

Introduced species (e.g., stages 2-5 after Colautti and 
Maclsaac 2004) are those transferred from their native 
range into a new nonindigenous area. Over the past cen¬ 
tury it has become increasingly clear how disruptive 
human-caused biological introductions have been to the 
planet. While not all introduced species cause obvious 
harm, some introduced species can eventually become 
economic threats and lead to serious conservation prob¬ 
lems (Simberloff et al. 1997). As of 2005, it was estimat¬ 
ed that the cost of enviromnental damages, losses, and 
control due to introduced species exceeded $120 billion 
per year in the United States alone (Pimentel et al. 2005). 
Prior to 2011, the state of Florida had 137 documented 
introduced reptile and amphibian taxa (56 being estab¬ 
lished), which ranks highest in the world (Krysko et al. 
2011a, 2012). Invasion pathways in Florida include (few¬ 


est to highest numbers) biological control, zoos, cargo/ 
plant shipments, and the pet trade. 

Pinesnakes, bullsnakes, and gophersnakes ( Pituophis 
Holbrook 1842) are large (up to 254 cm total length) 
constrictors native to North America, characterized by 
disproportionately small heads, four prefrontal scales, 
and a large rostral plate that extends upwards between 
the intemasals (Conant and Collins 1991). Based primar¬ 
ily on molecular data using Parsimony and Maximum 
Likelihood analyses with 893 base pairs (bp) of the nico¬ 
tinamide adenine dinucleotide dehydrogenase subunit 4 
(ND4) region (Rodriguez-Robles and De Jesus-Escobar 
2000), the P. melanoleucus species complex contains 
three currently recognized species; P. melanoleucus 
(sensu stricto; Pinesnakes; with three subspecies P. m. 
lodingi, P m. melanoleucus , P. m. mugitus ), P. cateni¬ 
fer (gophersnakes and bullsnakes; with six subspecies P. 
c. affinis, P. c. annectens, P c. catenifer , P. c. deserti- 
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Figure 1 . Map of Zoo Miami bounded in green. Note that major roadways, residential areas, and undeveloped protected lands 
surround zoo property. Dots represents locations of Pituophis found on zoo property; yellow = UF-Herpetology 157954 (gravid 
female) and red = UF-Herpetology 163092 (adult male). 


cola, P. c. pumilis, P. c. sayi ), and P ruthveni (Louisi¬ 
ana pinesnake). Pituophis melanoleucus (Daudin 1803) 
occurs in the eastern United States from southern New 
Jersey south to extreme southern peninsular Florida (i.e., 
Miami-Dade County; Krysko et al. 2011b) and west 
to Kentucky and southeastern Louisiana (Rodriguez- 
Robles and De Jesus-Escobar 2000). This species lacks 
a dark line from the eye to the angle of the jaw, has a 
dorsal pattern either absent (uniform black), obscure, or 
whitish to brownish with 23-30 distinct dark dorsal body 
blotches that are clearly separated from each other both 
anteriorly and posteriorly along the body and tail (Knight 
1986; Powell et al. 1998; Reichling 1995; Thomas et al. 
1976). Pituophis catenifer occurs from the Pacific Ocean 
east to Wisconsin, Illinois, and Texas, and from Canada 
south to Mexico (Rodriguez-Robles and De Jesus-Esco¬ 
bar 2000; Powell et al. 1998). This species typically has 
a dark line from the eye to the angle of the jaw, and a yel¬ 
low or cream-colored dorsal pattern with 41-79 distinct 
dark dorsal blotches that are clearly separated from each 
other both anteriorly and posteriorly along the body and 
tail (Knight 1986; Powell et al. 1998; Reichling 1995; 
Thomas et al. 1976). Pituophis ruthveni occurs in allo- 
patric populations in western-central Louisiana to eastern 
Texas (Ealy et al. 2004; Powell et al. 1998). This species 
sometimes lacks a dark line from the eye to the angle of 
the jaw, and has a pale brown dorsal pattern with 28-38 
dark dorsal blotches; near the head the blotches obscure 
the ground coloration, whereas near the tail they are dis¬ 
tinctly separated from each other (Knight 1986; Pow¬ 


ell et al. 1998; Reichling 1995; Stull 1929; Thomas et 
al. 1976). Although P. ruthveni is nested within a clade 
containing only P. c. sayi, it is recognized as a separate 
species because it occurs in allopatric populations and 
is somewhat diagnosable using color pattern characters 
(Collins 1991; Knight 1986; Reichling 1995; Rodriguez- 
Robles and De Jesus-Escobar 2000; Thomas et al. 1976). 
Pituophis ruthveni is also a candidate for listing as an 
imperiled species by the U.S. Fish and Wildlife Service 
(2013). 

The last known Pituophis melanoleucus from ex¬ 
treme southern peninsular Florida (UF-Herpetology 
45970) was collected in 1980 in a disturbed pineland 
(with Casuarina and Schinus) in Cutler Ridge, Miami- 
Dade County, and because of ongoing dense urbaniza¬ 
tion this species is believed to be extirpated along the 
Atlantic Coast Ridge (Krysko et al. 2011b). In 2010, two 
Pituophis were collected on the Atlantic Coast Ridge at 
Zoo Miami, Miami-Dade County; one was found in an 
undeveloped area and another near public access. Based 
on color pattern alone, these snakes were suspected to be 
introduced P. ruthveni and reported to represent the first 
known vouchers for this species in Florida (Krysko et al. 
2011a). Many documented introductions categorize spe¬ 
cies based on sometimes vague superficial morphology, 
such as color patterns, which may or may not be arbitrary 
human constructs. However, as more wide-range phylo¬ 
genetic studies are conducted and published, the oppor¬ 
tunity arises for other researchers to compare DNA from 
known data sets to suspected introduced individuals in 
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Verification of introduced Pituophis ruthveni in Florida 



Figure 2. Well-developed Pituophis embryo (UF-Herpetology 164295) oviposited from wild collected gravid female (UF-Herpe- 
tology 157954) in Miami, Miami-Dade County, Florida. 


order to confirm species identification as well as deter¬ 
mine their geographic origins. In this paper, we conduct 
molecular analyses of Pituophis in a coalescent frame¬ 
work to confirm species identity and phylogenetic place¬ 
ment of our two specimens, followed by more detailed 
examination of morphology and color pattern. 

Material and Methods 

Site description and specimen acquisition 

Zoo Miami is situated at 12400 SW 152 th Street, Miami, 
Miami-Dade County, Florida, USA (Fig. 1; 25.611926°N, 
80.398042°W, Datum WGS84, elev. 2 m). The property 
consists of ca. 300 ha, 106 ha of which are undeveloped 
managed lands, predominantly of pine rockland habitat. 
Zoo Miami property is surrounded by a mixture of natu¬ 
ral areas, disturbed areas, and a county park, followed by 
dense urbanization. 

On 16 May 2010 at 1645 h, an adult Pituophis (gravid 
female, 1,302 mm SVL, 1,486 mm TL; UF-Herpetolo¬ 
gy 157954; see Fig. 86 in Krysko et al. 2011) was col¬ 
lected in a service area behind a large animal exhibit 
(25.60395°N, 80.4006°W). This snake was observed 
by zoo staff the previous day along an adjacent public 
walkway, but was not captured. This snake was retained 
in captivity and oviposited three eggs on 22 June 2010. 
The eggs were viewed with a light on 28 June 2010; all 
three eggs contained a dark blood spot, but only one egg 
had an obvious network of veins developing. The first 
two eggs failed to develop and were discarded on 06 July 
2010. The third egg had an unpleasant odor and was fro¬ 
zen on 14 September 2010; it was dissected on 20 Sep¬ 


tember 2010 and revealed a well-developed embryo (UF- 
Herpetology 164295; Fig. 2). 

On 25 December 2010 at 1215 h, another adult (male, 
1,425 mm SVL, 1,635 mm TL) Pituophis (UF-Herpe¬ 
tology 163092) was collected in an undeveloped area 
(25.60304°N, 80.40295°W), across a large man-made 
lake and 0.26 km southwest of the first snake. 

The well-developed embryo, shed skins from the two 
adults, and digital images were deposited in the Divi¬ 
sion of Herpetology, Florida Museum of Natural His¬ 
tory, University of Florida. The female (UF-Herpetology 
157954) is currently housed at the Memphis Zoo, and the 
male (UF-Herpetology 163092) is housed at Zoo Miami. 

Laboratory techniques 

DNA isolations were obtained using QIAquick PCR Pu¬ 
rification Kit and DNeasy Blood and Tissue Kit (Qiagen 
Sciences, LLC). Using total cellular DNA as a template 
and Polymerase Chain Reaction (PCR) methodology 
(Saiki et al. 1988), mitochondrial DNA (mtDNA) was 
amplified and sequenced for the ND4 region follow¬ 
ing Rodriguez-Robles and De Jesus-Escobar (2000). 
The ND4 region includes a section of the 3 ’ end of the 
ND4 gene, and two subsequent transfer ribonucleic ac¬ 
ids (tRNA His , tRNA Ser ), which were sequenced using the 
primers ND4 and Leu (Arevalo et al. 1994). PCR was 
conducted in 25 pi reactions: 9.5 pi H 2 0, 12.5 pi GoTaq® 
Master Mix (Promega Corp, Madison, Wisconsin, USA), 
1.0 pi each primer (10 pM), and 1.0 pi DNA template. 
PCR parameters included initial denaturing at 94 °C for 
three min, followed by 35 cycles of amplification: de- 
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Figure 3. Maximum Likelihood phylogeny for Pituophis (Squamata: Colubridae) snakes, including the two known P. ruthveni 
(highlighted in yellow, UF-Herpetology 157954 and 163092) collected in Miami, Miami-Dade County, Florida. Note that values (> 
50%) above nodes represent bootstrap support. Inset photograph of UF-Herpetology 157954 by Dustin C. Smith. 


naturing at 94 °C for one min, annealing at 52 °C for 
one min, and extension at 72 °C for one min, followed 
by a final extension at 72 °C for seven min. Three pi of 
each PCR product were electrophoresed on a 1% agarose 
gel, visualized with GelRed™ staining (Biotium Inc., 
Hayward, California, USA), and compared with a DNA 
standard. Sequence files from the automated sequencer 
(Genomics Division, Interdisciplinary Center for Bio¬ 
technology Research, University of Florida) were as¬ 
sembled and edited as necessary with Geneious software 
(ver. 6.1, created by Biomatters. Available from http:// 
www.geneious .com). 

Phylogenetic analyses. —DNA sequence data were 
downloaded from GenBank for 46 snakes, including 42 
Pituophis , and one of each Lampropeltis getula, Panthe- 
rophis vulpinus , Bogertophis subocularis, and Arizona 
elegans incorporating the original data set from Rodri- 
guez-Robles and De Jesus-Escobar (2000) and current 
taxonomy after Pyron and Burbrink (2009). GenBank 
Accession numbers for our two Pituophis specimens 


(UF-Herpetology 157954 and 163092) are KJ938643 
and KJ938644, respectively. 

A total of 48 specimens with 875 base pairs (bp) of se¬ 
quence data were analyzed. Relationships among mtDNA 
haplotypes were estimated using both Maximum Likeli¬ 
hood (ML) and Bayesian Inference (BI) methods. ML 
was conducted with the General Time Reversible model 
with gamma distributed rate heterogeneity (GTR + T) 
and 1,000 nonparametric bootstrap replicates (Felsen- 
stein 2004) to assess node support using RAxML-HPC 
BlackBox (Stamatakis 2006; Stamatakis et al. 2008) 
from the CIPRES Science Gateway (Miller et al. 2010). 

BI was conducted using BEAST (ver. 1.8; Drummond 
and Rambaut 2007) from the UF-HPC Galaxy instance 
(http://hpc.ufl.edu; Blankenberg et al. 2010; Giardine et 
al. 2005; Goecks et al. 2010). The Bayesian Information 
Criterion in jModelTest (ver. 2.1.4; Darriba et al. 2012; 
Guindon and Gascuel 2003) determined the best-fit nu¬ 
cleotide substitution model to be Hasegawa, Kishino, 
and Yano with a proportion of invariant sites and gamma 
distributed rate heterogeneity (HKY + I + T). A relaxed 
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Figure 4. Bayesian Inference phylogeny for Pituophis (Squamata: Colubridae) snakes, including the two known P. ruthveni (high¬ 
lighted in yellow, UF-Herpetology 157954 and 163092) collected in Miami, Miami-Dade County, Florida. Note that values (> 95%) 
above nodes represent posterior probabilities. Inset photograph of UF-Herpetology 163092 by Dustin C. Smith. 


phylogenetics method was used without relying on a 
potentially arbitrary molecular clock (Zuckerkandl and 
Pauling 1965) that might incorporate uncertainty in the 
tree estimation process (Drummond et al. 2006). An un¬ 
correlated lognormal relaxed clock with coalescent con¬ 
stant population size (Kingman 1982), estimated base 
frequencies, randomly generated starting tree, and nor¬ 
mal distribution for the ucld.mean parameter priors were 
used. Two independent runs were performed consisting 
of three heated and one cold Markov Chain Monte Carlo 
(MCMC) estimated for 40 million generations, with ev¬ 
ery 1,000th sample being retained. Both MCMC runs 
were analyzed independently (to confirm chains were 
converging and not sampling local optima) using Tracer 
(ver. 1.6) for ESS values >200, as well as for a split stan¬ 
dard deviation less than 0.005 for -InL tree values among 
chains that indicate parameter stationarity was achieved. 
Trees sampled prior to stationarity were discarded as 
bum-in, which occurred prior to five million generations. 
Trees from both independent MCMC runs were com¬ 
bined and burn-in was removed using LogCombiner (ver. 
1.8), the best statistically supported tree (i.e., Maximum 
clade credibility tree) with mean heights was obtained 
using TreeAnnotator (ver. 1.8), and a phylogenetic hy¬ 


pothesis with posterior probabilities was created using 
FigTree (ver. 1.4). 

The most credible inferences of phylogenetic relation¬ 
ships were confined to nodes where nonparametric boot¬ 
strap values > 70% and posterior probability (Pp) was > 
95% (Hillis and Bull 1993; Felsenstein 2004). 

Morphology and color pattern 

We determined sex, snout-vent length (SVL), tail length, 
number of ventrals, subcaudals, supralabials, infralabi¬ 
als, preoculars, postoculars, temporals, loreals, and dor¬ 
sal scale rows; and color pattern of dorsum and venter. 
We compared these data to those found in the literature. 

Results 

Phylogenetic analyses. —Both ML and BI methods 
produced identical phylogenetic groupings (Figs. 3 and 
4). Although some of these clades are organized differ¬ 
ently in relation to one another the monophyly of Pituo¬ 
phis is well supported, which is congment with the find¬ 
ings by Pyron and Burbrink (2009), though the latter 
study used only single samples for each species. Both of 
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our two Pituophis specimens have the same mtDNA hap- 
lotype, and both phylogenetic methods place them within 
the P. catenifer sayi / P ruthveni clade. 

Morphological data for UF-Herpetology 157954 in¬ 
clude 226 ventrals, 55 subcaudals, 8/8 (left/right) supral- 
abials, 11/11 infralabials, 1/1 preoculars, 7/7 postoculars, 
4 temporals, 1/1 loreals, 27-30-24 dorsal scale rows, 
34 body blotches, 8 tail blotches, parietal stripe present, 
and heavily patterned venter. Data for UF-Herpetology 
163092 include 212 ventrals, 57 subcaudals, 8/8 (left/ 
right) supralabials, 11/11 infralabials, 1/1 preoculars, 7/7 
postoculars, 4 temporals, 1/1 loreals, 27-31-23 dorsal 
scale rows, 32 body blotches, and 11 tail blotches. 

Discussion 

Our ML and BI phylogenies produced identical mam 
phylogenetic groupings (Figs. 3 and 4) as those found 
in the ML analysis by Rodriguez-Robles and De Jesus- 
Escobar (2000). However, we found no support for some 
relationships, and no support values are provided on the 
original ML tree by Rodriguez-Robles and De Jesus- 
Escobar (2000). Our two Pituophis specimens were 
placed within a well-supported P. catenifer sayi / P ruth¬ 
veni clade, the same group of specimens (except for our 
Florida specimens) uncovered by Rodriguez-Robles and 
De Jesus-Escobar (2000). Pituophis catenifer sayi and 
P. ruthveni were also found to be sister taxa based on 
a combined mtDNA and single nuclear (nDNA) (Pyron 
and Burbrink 2009) and phenetic morphological similar¬ 
ity (Reichling 1995) analyses. Nonetheless, we found 
no molecular support for the recognition of P ruthveni 
as a separate species. One of the limitations of our and 
Rodriguez-Robles and De Jesus-Escobar’s (2000) mo¬ 
lecular studies is the use of only a single locus (ND4 re¬ 
gion), and adding additional mtDNA and unlinked nDNA 
genes might provide needed data for distinguishing be¬ 
tween these two named taxa. Pituophis ruthveni is cur¬ 
rently recognized as a separate species because it occurs 
in allopatric populations and is believed to be diagnos- 
able using color pattern characters, the most diagnostic 
being 28-38 dark dorsal body blotches and the blotches 
obscuring the ground coloration anteriorly (Collins 1991; 
Reichling 1995; Rodriguez-Robles and De Jesus-Esco- 
bar 2000). Our two Pituophis specimens exhibit these 
three characters, thus we categorized them as P ruthveni. 

Before our specimens were found, Pituophis ruthveni 
was not known to be kept at Zoo Miami, therefore this 
species is not representative of a zoo-mediated introduc¬ 
tion pathway and was likely released by an outside per¬ 
son. Other species such as the Reticulated python, Ma- 
layopython reticulatus (see Kaiser et al. 2013; Reynolds 
et al. 2014), and Pacific Coast giant musk turtle, Stauro- 
typus salvinii, are other examples of reptile species that 
have been illegally released on zoo property, the latter 
possibly established (Smith et al. 2011). Although we are 


currently uncertain if P. ruthveni is established in the vast 
protected undeveloped habitats surrounding public ac¬ 
cess areas, an adult male and gravid female were found 
suggesting reproduction might have taken place in the 
wild. 
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